Quantitative calculations of the excitonic energy spectra of semiconducting 
single-walled carbon nanotubes within a 7r-electron model 
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Using Coulomb correlation parameters appropriate for 7r-conjugated polymers (PCPs), and a 
nearest neighbor hopping integral that is arrived at by fitting the energy spectra of three zigzag 
semiconducting single- walled carbon nanotubes (S-SWCNTs), we are able to determine quantita- 
tively the exciton energies and exciton binding energies of 29 S-SWCNTs within a semiempirical 
TT-electron Hamiltonian that has been widely used for PCPs. Our work establishes the existence of 
a deep and fundamental relationship between PCPs and S-SWCNTs. 
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I. INTRODUCTION 

The photophysics of semiconducting single-walled 
carbon nanotubes (S-SWCNTs) are of strong cur- 
rent interest because of their potential applications in 
optoelectronics ii2, Recent theoretical works have empha- 
sized the Coulomb-induced binding between the optically 
excited electron and hole in these systems The 
excitonic energy spectra of S-SWCNTs are now under- 
stood qualitatively: there occur in these nanotubes a se- 
ries of energy manifolds labeled by an index n = 1, 2, 
... etc., with each manifold containing both optically al- 
lowed and dark exciton states, and a continuum band 
separated from the optical exciton by a characteristic ex- 
citon binding energy. Nonlinear optical absorption mea- 
surements have demonstrated that distinct energy gaps 
occur between the lowest two-photon allowed states and 
the n = 1 optical excitonsii2*ii*i2ii^ This energy gap is a 
lower bound to the binding energy of the n = 1 optical 
exciton^ Nonlinear absorption measurements involving 
the n = 1 continuum band^^ and the n — 2 exciton-'-^ 
have also been performed. 

Exciton formation in S-SWCNTs is a consequence of 
direct Coulomb electron-electron (e-e) interaction, taking 
proper account of which is a notoriously difficult many- 
body problem. Theoretical calculations of excitons and 
continuum band energies in the S-SWCNTs are therefore 
necessarily approximate. One semiquantitative approach 
that has been popular employs an ab initio approach for 
the ground state, which is followed by the determination 
of the quasiparticle energies within the GW approxima- 
tion and the solution of the Bethe-Salpeter equation of 
the two-particle Green's function,'''^ This approach be- 
comes difficult to implement for chiral S-SWCNTs with 
large diameters and large unit cells, precisely the sys- 
tems for which the nonlinear optical measurements are 
being performed ii2iiiii24iiiiii^ Scaling relationships for 
the exciton binding energies in wider chiral nanotubes 
have therefore been proposedAS Understanding nonlin- 
ear absorption within such indirect approaches is diffi- 
cult. 

An alternate approach to the exciton problem in 



the S-SWCNTs is to adopt a semiempirical 7r-electron 
approximation^^ that is widely used to describe planar 
TT-conjugated systems. The model assumes that the low- 
est energy excitations in planar conjugated systems in- 
volves the TT-electrons only and ignores the electrons oc- 
cupying orthogonal cr-bands. Since tt-tt* excitation en- 
ergies decrease rapidly with increasing system size, while 
the a and a* bands are nearly dispersionless, the ap- 
proximation is excellent for large systems. Thus tt- 
conjugated polymers (PCPs) have been discussed exten- 
sively within the semiempirical Pariser-Parr-Pople (PPP) 
Hamiltonian 4^ 
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Here c]^ creates a 7r-electron with spin a (t, i) in the Pz 
orbital of the ith carbon atom, (ij) implies nearest neigh- 
bor (n.n.) atoms i and j, ~ c[„Cic is the number of 
TT-electrons with spin a on the atom i, and Ui — is 
the total number of 7r-electrons on the atom. The param- 
eter tij is the one-electron hopping integral between Pz 
orbitals of n.n. carbon atoms, U is the on-site e-e repul- 
sion between two 7r-electrons occupying the same carbon 
atom Pz orbital, and is the intersite e-e interaction. 

The semiempirical model suffers from the disadvan- 
tage that determination of the parameters is commonly 
difficult. On the other hand, many-body problems that 
would be formidable within the ah initio approach, such 
as the enhancement of the ground state bond-alternation 
in polyacetylene by e-e interactions^ or the occurrence 
of the lowest two-photon state below the optical state in 
finite polyenes and polyacetylene fi^ can be understood 
within Eq. (1). The semiempirical approach should thus 
be taken as complementary to the ab initio one. 

The above observations were the basis of our calcula- 
tions within Eq. (1) of the electronic structures and exci- 
ton binding energies,^ and more recently, nonlinear opti- 
cal absorptions^^ in a limited number of S-SWCNTs, in 
spite of their nonplanarity. Justification of this procedure 
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can also be found in earlier work claiming that curva- 
ture effects of fuUerene molecules can be included within 
Eq. (1) by proper modifications of the parameters^ 
The c-e interaction and hopping integral parameters used 
in our calculations were the same ones that had been 
used earlier for the PCP poly-paraphenylenevinylene 
(PPV)iSi The calculated energy differences (~ 0.3-0.4 
eV) between dominant two-photon states observed in ul- 
trafast spectroscopy and the n = 1 optical excitons, in 
S-SWCNTs with diameters d > 0.8 nm, were remarkably 
close to what are observed experimentally.^^ Other exper- 
imental workiSiiiiiSii^ havc also determined the binding 
energy of the n = 1 optical exciton in S-SWCNTs with 
diameters 0.8-1.0 nm to be ^0.4 eV, in agreement with 
Ref. 5i Very recent work has also confirmed the predicted 
binding energy of the n = 2 exciton in two different S- 
SWCNTsi^ 



The calculated absolute energies of the n = I and 2 
excitons S-SWCNTs in our earlier work^^ii^ however, are 
much larger than the experimental values (by as much as 
0.5 eV). Equally importantly, comparisons of experimen- 
tal and calculated exciton binding energies were based on 
considerations of diameters alone: while the nanotubes 
probed experimentally are chiral with large diameters, 
the theoretical calculations were either for zigzag nan- 
otubes or for narrow chiral nanotubes, both with smaller 
unit cells. Direct comparisons of theory and experiment 
for the same systems were thus not possible. Finally, 
because the theoretical calculations could be performed 
for only a handful of nanotubes, family relationships, es- 
tablished experimentally,^^ could not be demonstrated 
within the theoretical model. 



In the present paper, we report calculated energy spec- 
tra of 29 S-SWCNTs, only 9 of which are zigzag, within 
Eq. (1). The diameters of the S-SWCNTs we consider 
range from 0.56 to 1.51 nm. These calculations have been 
possible because of substantive improvements in our com- 
putational techniques, while the improved results are a 
consequence of careful parametrization of the n.n. hop- 
ping integral Uj = t that takes into account curvature 
effects (see below). For each S-SWCNT we calculate the 
absolute energies En and E22 of the n — 1 and 2 op- 
tical excitons, and their binding energies E^i and Eb2, 
respectively. We compare all theoretical quantities to ex- 
perimentally determined ones i-^yi^^'^^i^^i^'^i^^'^^i^^i^? The 
large number of S-SWCNTs that could be considered al- 
lows us to investigate family relationships based on cal- 
culations of the energy ratio E22/E11 and to compare 
these data to experiments.^^ We find excellent agreement 
between the theory and experiments. Our work demon- 
strates convincingly that the photophysics of S-SWCNTs 
and PCPs can be understood within the same general 
theoretical framework, albeit with different hopping in- 
tegrals. 



TABLE I: Calculated and experimental (Ref. 22) n = 1 and 
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II. METHOD AND PARAMETRIZATION OF 
TT-ELECTRON MODEL 

We use the single configuration interaction (SCI) ap- 
proximation to compute the energies of the one electron- 
one hole excitations of the Hamiltonian in Eq. (1). The 
SCI approximation includes configuration mixings only 
between excitations that are singly excited from the 
Hartree-Fock (HF) ground state. The justifications for 
the approximation have been given before^*i^ 

We discuss next how the parametrizations within 
Eq. (1) were reached. Since e-e interactions depend 
only on the distance and not on the curvature of the 
S-SWCNTs, it is logical to parametrize the Vij exactly 
as in the PCPs^ viz., 

v., = , ^ (2) 
Ky'l + 0.6117i?2 

where i?y is the distance between carbon atoms i and j 
in A. The two free parameters in Eq. ^ are U and n; 
the latter is introduced to take into account the dielectric 
screening due to the mediumiSS U = 11.26 eV and k = 
1 correspond to the standard Ohno parametrization for 
finite molecules.^"' SCI calculations of the optical absorp- 
tion in PPV with five different U (between to 10 eV) 
and three different k (1, 2, and 3) had indicated that only 
with U = 8 eV and k = 2 was it possible to fit all four 
absorption bands at 2.4, 3.7, 4.7, and 6.0 eV in PPV. 
The same U and k were then used for quantitative calcu- 
lations of nonlinear and triplet absorptions in PPV. Ex- 
perimentally determined energies of the two-photon state 
that dominates nonlinear optical spectroscopjiS^ and the 
lowest triplet state2& agreed remarkably well with the 
theory. We have used the same U and k for S-SWCNTs. 

In the absence of any known procedure for determin- 
ing t in S-SWCNTs, we had earlieriSiiS, chosen the value 
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2.4 eV that is used for planar aromatic systemsiSi Curva- 
ture in the S-SWCNTs impHes smaUer tt-tt overlap be- 
tween n.n. C atoms and hence a smaller t. We arrive 
at the proper t by fitting the calculated En and £'22 
for three different zigzag nanotubes, (10,0), (13,0) and 
(17,0), against the corresponding experimental quanti- 
ties (see below for discussions of how the experimental 
quantities were arrived at by previous authors). The the- 
oretical exciton energies are for [/ = 8 eV and k = 2, and 
four difJerent t = 1.8, 1.9, 2.0, and 2.4 eV. Our resuhs 
are summarized in Table I, which clearly indicates that 
t — 2.4 eV is too large and considerably better fits are 
obtained with t = 1.8-2.0 eV. We have chosen t — 2.0 eV 
for the complete set of calculations we report below. The 
fits in Table I improve with increasing nanotube diame- 
ter, implying that strictly speaking the hopping integral 
is diameter-dependent even within the range of diame- 
ters we consider. We do not attempt further fine tuning 
of parameters as this would necessarily lead to loss of 
simplicity and generality. 

As in previous workj^*i^ we use open boundary con- 
dition (OBC). Surface states due to danghng bonds at 
the nanotube ends appear in the HF band structure, 
and are discarded at the SCI stage of our calculations. 
The chiral S-SWCNTs we investigate have gigantic unit 
cells. The number of unit cells we retain depend both 
on the size of the unit cell and the convergence behavior 
of Ell. The procedure involved calculating the standard 
tight-binding (TB) band-structure with periodic bound- 
ary condition (PEC), and then comparing the PBC En 
with that obtained using OBC with a small number of 
unit cells. The number of unit cells in the OBC calcula- 
tion is now progressively increased until the difference in 
the computed En between OBC and PBC is less than 
0.004 eV (worst case). It is with this system size that the 
SCI calculations are now performed using OBC. Thus for 
example, our calculations for (7,0), (6,4), and (7,5) SWC- 
NTs are for 70, 16, and 5 unit cells, respectively, con- 
taining 1960, 2432, and 2180 carbon atoms, respectively. 
Since energy convergences are faster in the calculations 
with nonzero e-e interactions than the calculations in the 
TB limit, we are confident that this procedure gives ac- 
curate results. We retain an active space of 100 valence 
and conduction band states each in the SCI calculations. 
Stringent convergence tests involving gradual increase in 
the size of the active space indicate that the computa- 
tional errors due to the energy cutoff is less than 0.005 
eV (worst case). In addition to En and E22 we also 
calculate the corresponding exciton binding energies Ei,i 
and Eii2 , with the n ~ 1 and 2 continuum band threshold 
energies defined within the SCI to be the corresponding 
HF gaps. 



III. RESULTS AND DISCUSSION 

In Table II we have listed our calculated En, E22, 
Ebi, and Eb2 for 29 S-SWCNTs. We compare each 



of these quantities to those obtained by experimental 
investigators Nearly half the exciton en- 
ergies listed in Table II as experimental were obtained 
directly from spectroflurometric measurementaSi or from 
resonant Raman spectroscopy.^^ Using the experimental 
data in Ref. |^ Weisman and BachiloS^ derived empir- 
ical equations for the exciton energies of nanotubes for 
which direct experimental information do not exist. The 
experimental En and £^22 in Table I are obtained from 
these empirical equations. Dukovic et alm^ have given an 
empirical equation for the binding energy of the n = 1 
exciton, which was also derived by fitting the set of Ei,i 
obtained from direct measurements. Only two measured 
values exist currently for the binding energy of n = 2 
excitoui^ We make distinctions between the experimen- 
tal and empirical data in Table II, but in the text below 
we refer to both as experimental quantities. In Fig. 1 we 
have plotted the theoretical and experimental En and 
E22 against 1/d, while in Fig. 1 inset we show the errors 
in our calculations, AEn and AE22, defined as the cal- 
culated energies minus the experimental quantities. As 
seen in the figure the theoretical plots are nearly inde- 
pendent of chirality, and depend primarily on diameter. 
This is a consequence of Eq. (1), within which the ener- 
getics depend only on the conductivity. We find excellent 
agreement between calculated and experimental £11 for 
d > 0.75 nm, with |Ai?ii| < 0.1 eV. (The black arrow on 
the X-axis in Fig. 1 indicates d = 0.75 nm.) The agree- 
ment for d > 1 nm is even better with |Ai?ii| < 0.05 
eV. The larger disagreement with experiment (and the 
greater chirality dependency of the experimental ener- 
gies) for d < 0.75 nm is due to the breakdown of the tt- 
electron approximation. The disagreements between cal- 
culated and experimental E22 are larger, but even here 
the magnitude of the maximum error for d > 0.75 nm 
is within 0.2 eV, which is the C-C bond stretching fre- 
quency that can influence experimental estimation of ex- 
citon energies.^ The origin of the larger disagreement in 
the n — 2 region is the SCI approximation and not the 
TT-electron model: inclusion of higher order CI is more 
important in general for higher energy states. 

The agreement between the calculated and the experi- 
mental exciton binding energies in Table II is even more 
striking than the fits to the absolute energies. The dis- 
crepancies between theory and experiment is less than 
10%, which is the uncertainty in the empirical binding 
energiesiii Both E^i and Eb2 are inversely proportional 
to the diameter d and can be fitted approximately by 

Ebi ~ ^ eV, (3a) 

Eb2^ ^eV. (3b) 

The exciton binding energies depend weakly on chiral- 
ity. For both binding energies, the fits of Eq. (3) become 
better for larger diameter nanotubes. Eq. (3a) is remark- 
ably close to the empirical formula E^i ~ eV given 
in previous experimental workiii Although there are not 
enough experimental data to verify the £{,2 relation in 
Eq. (3b), the linear dependence against 1/d should be 
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TABLE II: Comparison of calculated and experimental/empirical n = 1 and 2 exciton energies and binding energies. The 
empirical exciton energies (Ref. 22) and exciton binding energies (Ref. El are in parentheses. 
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true for both binding energies. The calculated Ebi and 
i?fc2 are very close to those obtained earlier by us for the 
wide nanotubes ■with t = 2.4 eVf^*i^ even as the calcu- 
lated Ell and E22 are now quite different. This is simply 
a consequence of localization of the electrons at the large 
U /t considered here: exciton binding energies in this case 
depend primarily on electron correlation. Recall, for ex- 
ample, that in the limit where the intersite Coulomb in- 
teraction is cut off beyond the n.n. interaction Vi, the 
exciton binding energy is determined almost entirely by 
the difference in Coulomb interactions [/ — Vi,^^ with t 
playing a weak role for U /t > 4. 

Our ability to calculate En and E22 for a large num- 
ber of nanotubes allows us to demonstrate family be- 



havior within Eq. (1), which was not possible before. It 
has been shown that in plots of experimental E22/E11 
against the wavelength A22 corresponding to i?22, fami- 
lies of (n, m) S-SWCNTs with the same n — to lie on the 
same continuous curves^ The experimental energy ra- 
tios diverge from a central limiting number ~1.75, with 
the deviations increasing with n — m, and with the di- 
rection of the deviation being positive for the families 
(n — to) mod 3=1 and negative for the families (n — to) 
mod 3 = 2, respectively!^ In Fig. 2 we have plotted our 
calculated E22/E11 against the calculated A22- The cal- 
culated ratios do not show the same smooth behavior as 
in Ref. 21, since as already noticed in Fig. 1 inset our 
errors with E22 are larger than those with En. Never- 
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FIG. 1: (Color online) Calculated (red solid line, circle and square symbols) vs. experimental (blue dotted line, diamond and 
triangle symbols) En and £22 for 29 S-SWCNTs. The inset shows errors AEu {i = 1, 2) in the calculations, defined as the 
calculated minus the experimental or empirical energies. 
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FIG. 2: (Color online) Calculated E22/E11 vs. the wave- 
length A22 corresponding to £22. The blue solid line with 
circle symbols (red dashed line with triangle symbols) denote 
families (n — m) mod 3 = 1 (2). 



theless, two different classes of behavior, for the families 
(n — to) mod 3 = 1 and {n — to) mod 3 = 2 are very 
clear. The difference between the largest and the smallest 
ratios in the narrow diameter (small A22) end is smaller 
than what is seen experimentally, but this is once again 
merely a consequence of the larger errors in our calcula- 
tions for the narrower nanotubes. More importantly, the 



limiting E22/E11 in the wide diameter (large A22) end in 
Fig. 2 is very close to the experimental limiting ratio of 
1.75. 



IV. CONCLUSION 

In summary, direct calculations of the energy spectra 
of a large number of S-SWCNTs within the 7r-electron 
PPP model Hamiltonian and with a fixed set of param- 
eters show excellent agreement with experiments. The 
magnitudes of two of the three free parameters, U and 
K for Coulomb interactions, are the same as those used 
earlier to calculate the energy spectrum of PPVj^S the 
third free parameter t is obtained by fitting the energet- 
ics of three zigzag S-SWCNTs to account for the curva- 
ture effects. The calculated exciton energies and exciton 
binding energies in the first two manifolds are in excellent 
agreement with experiments, and the "family behavior" 
is demonstrated. In the diameter range of our calculated 
S-SWCNTs, we found linear dependence against \/d for 
the binding energies of excitons in both manifolds. It 
is unlikely that the agreement between theory and ex- 
periments is fortuitous. The Coulomb parameters used 
here were obtained after extensive search through a very 
large parameter space for PPV, and the same parameters 
successfully reproduced the energies of the dominant two- 
photon state^^ and the triplet energy spectrum'^^ of PPV. 
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Similarly, slightly smaller hopping in S-SWCNTs than in 
PPV is to be expectediiSi Our work demonstrates a uni- 
versahty in the photophysics of S-SWCNTs and PCPs 
that arises from their common quasi-one-dimensionality 
and TT-conjugation. The existing rich literature on the 
photo- and device physics of PCPs can therefore provide 
valuable guidance in the search for optoelectronic appli- 



cations of S-SWCNTs. 
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